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Abstract 

The nolinear hydro-dynamic equations of the surface of a liquid drop are 
shown to be directly connected to Korteweg de Vries (KdV, MKdV) systems, 
giving traveling solutions that are cnoidal waves. They generate multiscale 
patterns ranging from small harmonic oscillations (linearized model), to non- 
linear oscillations, up through solitary waves. These non-axis-symmetric lo- 
calized shapes are also described by a KdV Hamiltonian system. Recently 
such "rotons" were observed experimentally when the shape oscillations of a 
droplet became nonlinear. The results apply to drop-like systems from cluster 

formation to stellar models, including hyperdeformed nuclei and fission. 
47.55.Dz, 24.10.Nz, 36.40.-c, 97.60.Jd 



1 



Typeset using REVTfjX 



A fundamental understanding of non-linear oscillations of a liquid drop (NLD), which 
reveals new phenomena and flows more complicated than linear theory suggests, is needed in 
diverse areas of science and technology. Besides their direct use in rheological and surfactant 
theory such models apply to cluster physics 0, super- and hyper-deformed nuclei [|TJ, 

nuclear break-up and fission thin films ||, radar and even stellar masses and 

supernova [|T1,|To| . Theoretical approaches are usually based on numerical calculations within 
different NLD models, |2|-§J and explain/predict axis- symmetric, non-linear oscillations that 
are in very good agreement with experiment [l|J^-0|. However, there are experimental results 
which show non-axis-symmetric modes; for example, traveling rotational shapes |5]|| that 
can lead to fission, cluster emission, or fusion |5] [7J . 

In this letter the existence of analytic solutions of NLD models that give rise to traveling 
solutions which are solitary waves is proven. Higher order non-linear terms in the devi- 



ation of the shape from a sphere produce surface oscillations that are cnoidal waves |TT 
By increasing the amplitude of these oscillations, the non-linear contribution grows and 
the drop's surface, under special conditions (non-zero angular momentum), can transform 
from a cnodial wave form into a solitary wave. This same evolution can occur if there is 
a non-linear coupling between the normal modes. Thus this approach leads to a unifying 
dynamical picture of such modes; specifically, the cnoidal solution simulates harmonic os- 
cillations developing into anharmonic ones, and under special circumstances these cnoidal 
wave forms develop into solitary waves. Of course, in the linear limit the theory reproduces 
the normal modes of oscillation of a surface. 

Two approaches are used: Euler equations fHH, and Hamiltonian equations, which de- 
scribe the total energy of the system 0. We investigate finite amplitude waves, for which 
the relative amplitude is smaller than the angular half-width. These excitations are also 
"long" waves, important in the cases of externally driven systems, where the excited wave- 
length depends by the driving frequency. The first original observations of travelling waves 
on liquid drops are described in [Q . Similar travelling or running waves are also discussed or 
quoted in |2|||. These results suggest that higher amplitude non-linear oscillations can lead 



to a traveling wave that originates on the drop's surface and developes towards the interior. 
This is shown to be related in a simply way to special solitary wave solutions, called "rotons" 



in the present analysis. Recent experiments and numerical tests PJlzfl suggest the existence 
of stable traveling waves for a non-linear dynamics in a circular geometry, re-enforcing the 
theory. 

A new NLD model for describing an ideal, incompressible fluid drop exercising irrota- 
tional flow with surface tension, is employed in the analysis. Series expansion in terms of 
spherical harmonics are replaced by localized, nonlinear shapes shown to be analytic solu- 
tions of the system. The flow is potential and therefore governed by Laplace's equation for 
the potential flow, A$ = 0, while the dynamics is described by Euler's equation, 

p(d t v + (v-V)v) = -VP + f } (1) 

— * 

where P is pressure. If the density of the external force field is also potential, / = — V\l/ 
where \I/ is proportional to the potential (gravitational, electrostatic, etc.), then Eq. ([[]) 
reduces to Bernoulli's scalar equation. The boundary conditions (BC) on the external free 
surface of the drop, SI, and on the inner surface S2, [P,P|JfT|, are f| si = (r t + r e 9 + 
an d ^|s2 = 0, respectively. <3> r = f is the radial velocity, $g = r 2 9, $^ = r 2 sin 00 
are the tangential velocities. The second BC occurs only in the case of fluid shells or 
bubbles. A convenient geometry places the origin at the center-of-mass of the distributon 
r{9, 0, t) = Rq[1 + g(9)rj((j) — Vt)] and introduces for the dimensionless shape function grj a 
variable denoted £. Here Rq is the radius of the undeformed spherical drop and V is the 
tangential velocity of the traveling solution £ moving in the direction and having a constant 
transversal profile g in the 9 direction. The linearized form of the first BC, f|si = 7"t|si, 
allows only radial vibrations and no tangential motion of the fluid on SI, |2|,|3|,[11[] . The second 
BC restricts the radial flow to a spherical layer of depth h{9) by requiring $ r \r=Ro-h = 0- 
This condition stratifies the flow in the surface layer, Rq — h < r < Ro(l + £), and the 
liquid bulk r < R — h. In what follows the flow in the bulk will be considered negligible 
compared to the flow in the surface layer. This condition does not restrict the generality of 



the argument because h can always be taken to be Ro. Nonetheless, keeping h < Ro opens 
possibilities for the investigation of more complex fluids, e.g. superfluids, flow over a rigid 
core, multilayer systems |||7|] or multiphases, etc. Instead of an expansion of $ in term of 
spherical harmonics, consider the following form 

$ = £(r/itW) n /n(0,<M)- (2) 

n=0 

The convergence of the series is controlled by the value of the small quantity e = max\ |, 
|TT| . The condition max\h/Ro\ ~ e is also assumed to hold in the following development. 
Laplace's equation introduces a system of recursion relations for the functions f n , = 
—fi — Aq/ /2, etc., where Aq is the (9, </>) part of the Laplacean. Hence the set of unknown 
/ n 's reduces to f and fi. The second BC, plus the condition £^ = — V£ t , for traveling waves, 
yields to second order in e, 

/ M = VRlsm 2 9£(l + 20 /h + O 3 (0, (3) 

i.e., a connection between the flow potential and the shape, which is typical of nonlinear 
systems. Eq.(3) together with the relations fx ~ R^ t ~ — — ' w hich follow from 
the BC and recursion, characterize the flow as a function of the surface geometry. The 
balance of the dynamic and capillary pressure across the surface SI follows by expanding 
up to third order in £ the square root of the surface energy of the drop 

Us = aRl Jjl + 0^/(1 + 2 + # + sin 2 (4) 

and by equating its first variation with the local mean curvature of SI under the restriction 
of the volume conservation. The surface pressure, in third order, reads 

P|ei = |-(-2£ - 4£ 2 - Api + 3te e ctg9), (5) 

where a is the surface pressure coefficient and the terms and £0,0 are neglected 

because the relative amplitude of the deformation e is smaller than the angular half-width 
L, £ = i<t><f> — t 2 /L 2 <^ 1, as most of the experiments P,|T|,|9|,|T2|] concerning traveling surface 



patterns show. Eq.(5) plus the BC yield, to second order in e, 



V 2 RUm 2 9~ 
$ U. _| o £ 2 



rr (2£ + 4£ 2 + A n £-3£V^). (6) 



pi? 

The linearized version of Eq. (||) together with the linearized BC, $ r |si = -Ro£t, yield a 
limiting case of the model, namely, the normal modes of oscillation of a liquid drop with 
spherical harmonic solutions 0,0]. Differentiation of Eq. (H) with respect to <p together with 
Eqs.(3,5) yields the dynamical equation for the evolution of the shape function r/(0 — Vt): 

Arjt + Br]^ + Cgr]r]^ + Drj^ = 0, (7) 



which is the Korteweg-de Vries (KdV) equation [] L] with coefficients depending parametri- 
cally on 9 

m{Ro + 2h} sin 2 9 a (2g + A^) 
A=V , B = , 

h pRo g 

C = 8 C R 8h^ ~~^)' D = ~pR lin 2 9- (8) 

In the case of a two-dimensional liquid drop, the coefficients in Eq. (|8]) are all constant. 
Eq. (0) has traveling wave solutions in the direction if Cg/(B — AV) and D/(B — AV) 
do not depend on 9. These two conditions introduce two differential equations for g{9) 
and h(9) which can be solved with the boundary conditions g = h = for 9 = 0, it. For 
example, hi = Rq sin 2 9 and g\ = P 2 {9) is a particular solution which is valid for h Rq. It 
represents a soliton with a quadrupole transvere profile, being in good agreement with 
The next higher order term in Eq. (H), —3^ 2 ^ectg9, introduces a rfr]^ nonlinear term into 



the dynamics and transforms the KdV equation into the modified KdV equation [11]. The 



traveling wave solutions of Eq. (0) are then described by the Jacobi elliptic function ITT 



V = a 3 + (a 2 - a 3 )sn 2 (y [ 2D (<f> ~ Vt); m) , (9) 

where the aij are the constants of integration introduced through Eq. (|7|) and are related 
through the velocity V = C{a x + a 2 + a 3 )/3A + B/A and m 2 = gE^. m G [0, 1] is the 
free parameter of the elliptic sn function. This result for Eq. (PI) is known as a cnoidal wave 



solution with angular period T = K[m] yC(«3 — a\)/3D where K{m) is the Jacobi elliptic 
integral. If a 2 — > oci — > 0, then m — ► 1, T — > oo and a one-parameter (770) family of traveling 
pulses (solitons or anti-solitons) is obtained, 

ilsoi = ilosech 2 [{<p-Vt)/Ll (10) 

with velocity V = r] C/3A + B/A and angular half- width L = ^jviD jCr\Q. Taking for the 
coefficients A to D the values given in Eq.(8) for 9 = tt/2 (the equatorial cross section) 
and hi, gi from above, one can calculate numerical values of the parameters of any roton 
excitation function of r)o only. 

The soliton, among other wave patterns, has a special shape-kinematic dependence r] ~ 
V ~ 1/L; a higher soliton is narrower and travels faster. This relation can be used to 
experimentally distinguish solitons from other modes or turbulence. When a layer thins 
(h — > 0) the coefficient C in eq.(8) approaches zero on average, producing a break in the 
traveling wave solution (L becomes singular) because of the change of sign under the square 
root, eq.(9). Such wave turbulence from capillary waves on thin shells was first observed 
in [9]. For the water shells described there, eq.(8) gives h(fim) < 20u/k, that is h=15- 
25yum at V=2.1-2.5ms _1 for the onset of wave turbulence, in good agreement with the 
abrupt transition experimentaly noticed. The cnoidal solutions provide the nonlinear wave 
interaction and the transition from competing linear wave modes (C < 0) to turbulence 
(C ~ 0). In the KdV eq.(7), the nonlinear interaction balances or even dominates the linear 
damping and the cnoidal (roton) mode occurs as a bend mode (h small and coherent traveling 
profile) in agreement with [9] . The condition for the existence of a positive amplitude soliton 
is gCD > which, for g < 0, limits the velocity from below to the value V > huj2/Ro where 
uj 2 is the Lamb frequency for the A = 2 linear mode @||. This inequality can be related 
to the "independent running wave" described in 0, which lies close to the A = 2 mode. 
Moreover, since the angular group velocity of the (A,/i) normal mode, V\ tfl = tO\j fi, has 
practically the same value for A = 2 (yU = 0,±1, tesseral harmonics) and for A = /x, any 
A (sectorial harmonics) this inequality seems to be essential for any combination of rank 2 
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tesseral or sectorial harmonics, in good agreement with the conclusions in . The periodic 
limit of the cnoidal wave is reached for m ~ 0, that is, 02 — 03 — 0, and the shape is 
characterized by harmonic oscillations (sn —>■ sin in Eq. (0)) which realize the quadrupole 
mode of a linear theory Y2 limit [|2],|3|] or the oscillations of tesseral harmonics ||, Fig. 1. 

The NLD model introduced in this paper yields a smooth transition from linear oscilla- 
tions to solitary traveling solutions ("rotons") as a function of the parameters ati; namely, 
a transition from periodic to non-periodic shape oscillations. In between these limits the 
surface is described by nonlinear cnoidal waves. In Fig.l the transition from a periodic limit 
to a solitary wave is shown, in comparison with the corresponding normal modes which can 
initiate such cnoidal nonlinear behavior. This situation is similar to the transformation of 
the flow field from periodic modes at small amplitude to traveling waves at larger amplitude 
PIJ. The solution goes into a final form if the volume conservation restriction is enforced: 
/ E (l + g(9)rj((j),t)) 3 dfl = An and requires T](4>,t) to be periodic. The periodicity condition, 
nK\{aLs — 012) / i a 3 ~ a i)] = fty/o^ ~ a i f° r an Y positive integer n, is only fulfilled for a finite 
number of n values, and hence a finite number of coresponding cnoidal modes. In the roton 
limit the periodicity condition becomes a quasi-periodic one because the amplitude decays 
rapidly. This approach could be extended to describe elastic modes of surface as well as 
their nonlinear coupling to capillary waves. The double-periodic structure of the elliptic 
solutions could describe the new family of normal wave modes predicted in [[| . 

The development up to this point was based on Euler's equation. The same result 
will now be shown to emerge from a Hamiltonian analysis of the NLD system. Recently, 
Natarajan and Brown 0] showed that the NLD is a Lagrangian system with the volume con- 
servation condition being a Lagrange multiplier. In the third order deviation from spherical, 
the NLD becomes a KdV infinite-dimensional Hamiltonian system described by a nonlin- 
ear Hamiltonian function H = J 27r H,d(f). In the linear approximation, the NLD is a linear 
wave Hamiltonian system |2|rj] . If terms depending on 9 are absorbed into definite integrals 
(becoming parameters) the total energy is a function of 77 only. Taking the kinetic energy 
from PHI, $ from Eq. (g) and using the BC, the dependence of the kinetic energy on the 
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tangential velocity along 9 direction, becomes negligible and the kinetic energy can be 
expressed as a T[rj\ functional. For traveling wave solutions d t = —Vd^, to third order in e, 
after a tedious but feasible calculus, the total energy is: 

E = J*\c lV + C 2V 2 + C 3V 3 + C^l)d<f>, (11) 

where C x = 2afi*Sl$, C 2 = aR 2 ^ + S$/2) + R 6 pV 2 C 3 2 ;Zl/2, C 3 = aR 2 S^/2 + 
R 6 pV 2 (2S 3 Jl 1 2 R + S 5 !^ + RoS*l£)/2, C 4 = <rB%S$°/2, with Sfj = Rq 1 f£ h l g i g J e sia k edd. 
Terms proportional to rji] 2 can be neglected since they introduce a factor tjq/L 2 which is 
small compared to t]q, i.e. it is in the third order in e. If Eq. ([□]) is taken to be a Hamilto- 
nian, E — > H[rj], then the Hamilton equation for the dynamical variable rj, taking the usual 
form of the Poisson bracket, gives 

rjtdcj) = (2C 2 r]<t> + 6C3VV<f> ~ 2C^i] H<i ) d<f>. (12) 

Since for the function r)(<fi — Vt) the LHS of Eq.(12) is zero, the integrand in the RHS gives 
the KdV equation. Hence, the energy of the NLD model, in the third order, is interpreted as 
a Hamiltonian of the KdV equation f?,II|. This is in full agreement with the result finalized 
by Eq. (0) for an appropriate choice of the parameters and the Cauchy conditions for g, h. 
The dependence of E(a\, a 2 )\voi=constant, Eq.(ll), shows an energy minimum in which the 



solitary waves are stable, [12 



The nonlinear coupling of modes in the cnoidal solution could explain the occurence of 
many resonances for the I = 2 mode of rotating liquid drops, at a given (higher) angular 
velocity, fL3f . The rotating quadrupole shape is close to the soliton limit of the cnoidal wave. 
On one hand, the existence of many resonances is a consequence of by the multi- valley profile 
of the effective potential energy for the KdV, (MKdV) equation: rfi. — ar] + brf + erf + (rfr^ 4 ), 
fTT| . The frequency shift predicted by Busse in [[HJ can be reproduced in the present theory 



by choosing the solution h\ = Rosin9/2. It results the same additional pressure drop in the 
form of V 2 pR^sin 2 9 /2 like in [IB], and hence a similar result. For a roton emerged from a I = 
2 mode, by calculating the half-width (L 2 ) and amplitude (r) ma x j2 ) which fitt the quadrupole 
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shape it results a law for the frequency shift: Au 2 /u)2 — (1 ± 4L 2 («3 — a2)/3Ro)~ 1 V/u>2, 
showing a good agreement with the observations of Annamalai et al in i.e. many 
resonances and nonlinear dependence of the shift on = V. The special damping of the 
I = 2 mode for rotating drops could also be a consequence of the existence of the cnoidal 
solution. An increasing in the velocity V produces a modification of the balance of the 
coefficients C/D which is equivalent with an increasing in dispersion. 

The model introduced in this article proves that traveling analytic solutions exist as 
cnoidal waves on the surface of a liquid drop. These traveling deformations ("rotons") can 
range from small oscillations (normal modes), to cnoidal oscillations, and on out to solitary 
waves. The same approach can be applied to bubbles as well, except that the boundary 
condition on E 2 is replaced by a far-field condition [2,3] (recently important in the context 
of single bubble sonoluminiscence) . Nonlinear phenomena can not be fully investigated with 
normal linear tools, e.g. spherical harmonics. Using analytic non-linear solutions sacrifices 
the linearity of the space but replaces it with multiscale dynamical behavior, typical for non- 
linear systems (solitons, wavelets, compactons []T2|). They can be applied to phenomena like 
cluster formation in nuclei, fragmentation or cold fission, the dynamics of the pellet surface 
in inertial fusion, stellar models, and so forth. 

Supported by the U.S. National Science Foundation through a regular grant, No. 
9603006, and a Cooperative Agreement, No. EPS-9550481, that includes matching from 
the Louisiana Board of Regents Support Fund. 



9 



REFERENCES 
[1] R. E. Apfel et al, Phys. Rev. Lett. 78, 1912 (1997). 

[2] R. Natarajan and R. A. Brown, J. Fluid Mech. 183, 95 (1987); Phys. Fluids 29, 2788 
(1986); J. A. Tsamopoulos and R. A. Brown, J. Fluid Mech. 127, 519 (1983). 

[3] H-L. Lu and R. E. Apfel, J. Fluid Mech., 222, 351 (1991); T. Shi and R. E. Apfel, Phys. 
Fluids 7, 1545 (1995). 

[4] Y. Tian, R. G. Holt and R. E. Apfel, Phys. Fluids 7, 2938 (1995). 

[5] E. H. Trinh, P. L. Marston and J. L. Robey, J. Colloid Interface Sci. 124, 95 (1988). 

[6] E. Trinh and T. G. Wang, J. Fluid Mech. 122, 315 (1982). 

[7] E. H. Trinh, R. G. Holt, and D. B. Thiessen, Phys. Fluids 8, 43 (1995); P. L. Marston 
and S. G. Goosby, Phys. Fluids 28, 1233 (1985). 

[8] A. Ludu et al, Int. J. Modern Phys. E 2, 855 (1993); J. Phys. G: Nucl. Part. Phys. 21, 
1715 (1995). 

[9] R. G. Holt and E. H. Trinh, Phys. Rev. Lett. 77, 1274 (1996). 
[10] R. H. Durisen et al Astrophys. J. 305, 281 (1996). 

[11] G. L. Lamb, Elements of Soliton Theory, (John Wiley & Sons, New York, 1980). 

[12] Y. G. Kevrekidis; A. Ustinov; A. Ludu and J. P. Draayer, TheCNLS 17 th Ann. Conf, 
May 12-16, 1997, Los Alamos; P. Rosenau and J. M. Hyman, Phys. Rev. Let. 70, 564 
(1993); S. Dramanyan, et al, Phys. Rev. E 55, 7662 (1997). 

[13] F. H. Busse, J. Fluid. Mech. 142, 1 (1984); P. Annamalai, E. Trinh and T. G. Wang, 
J. Fluid. Mech. 158, 317 (1985). 



10 



FIGURES 

FIG. 1. The cnoidal solution for = 0. The soliton limit and a 3- and 4-mode solution is 
shown. The closest spherical harmonics to each of the cnoidal wave profiles (labelled Cn and Sol, 
respectively) is given for comparison. The labels A, [i and the parameters ai,2,3 of the coresponding 
cnoidal solution are given. 
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